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Abstract 

CN \ The couphngs of pions with lowest three doublets (0^, 1^), (0+, 1+) and (l+,2+) 

■ of heavy mesons are studied with light-cone QCD sum rules in the leading order of 

^ . heavy quark effective theory. The ambiguity due to presence of two distinct 1"*" states 

(N ■ are solved. 
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"Sh! 1 Introduction 

ix. 

^ • Remarkbale progress has been made in understanding the physics of heavy mesons com- 

posed of a heavy quark and a light quark with the discovery of the heavy quark symmetry. 
To obtain detailed predictions one has to employ some specific nonperturbative methods, 
^ I among which QCD sum rules [|l| is very useful. The heavy quark effective theory (HQET) 

' 1^ provides a systematic expansion of QCD in terms of l/mq, where mg is the heavy 

quark mass. The spectrum of the ground state heavy meson has been studied with the 
QCD sum rules in HQET in In the mass of the lowest excited heavy meson doublets 
(2+, 1'*') and (l"*", 0+) were studied with QCD sum rules in the heavy quark effective theory 
(HQET) up to the order of C(l/mQ). 

QCD sum rules has been used to analyse the couplings of the heavy hadrons with 
pions [5-15]. The widths for pionic decays of the lowest two excited doublets (2+, 1+) and 
(1"'",0"'") is calculated with conventional QCD sum rules in the leading order of HQET in 
0. Light-cone QCD sum rules (LCQSR) with the full QCD Lagrangian is first employed 
to derive the strong coupling constants: gnD'D, QttB'B [HIl- Recently the strong coupling 



constants: g^B^B*^ 9-nB[B[^ and g^B'^B* were calculated using LCQSR with finite heavy quark 
mass [|l3l. The couplings of heavy baryons with soft pions have been estimated with QCD 
sum rules in an external axial field |T^. In this work we employ the LCQSR in HQET to 
calculate the on-shell couplings of the pions with heavy meson doublets (0~, 1~), (0^, l'^) 
and (1+, 2+) to the leading order of 

The LCQSR is quite different from the conventional QCD sum rules, which is based 
on the short-distance operator product expansion (OPE). The LCQSR is based on the 
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OPE on the light cone, which is the expansion over the twists of the operators. The 
main contribution comes from the lowest twist operator. Matrix elements of nonlocal 
operators sandwiched between a hadronic state and the vacuum defines the hadron wave 
functions. When the LCQSR is used to calculate the coupling constant, the double Borel 
transformation is always invoked so that the excited states and the continuum contribution 
can be treated quite nicely. Moreover, the final sum rule depends only on the value of the 
wave function at a specific point like fniuo = which is much better known than the 
whole wave function |]Tl] . 



One difficult problem encountered in studying the decay widths of excited heavy mesons 
with QCD sum rules is the following. Except for the lowest states 0~, 1~, the spectra 
contains a pair of states for any spin-parity with close values in their masses but quite 
different in magnitudes of their decay widths. In our case, one of the two l"*" states is a 
narrow resonance decaying mainly by emitting a D wave pion, while the other one is a very 
wide resonance decaying by emitting a S wave pion. An interpolating current used for the 
narrow l"*" state with a small coupling to the other 1+ state may cause sizable error in the 
result of calculation. It is only in the rriQ oo limit, there is a conserved quantum number 
je, the angular momentum of the light component, which can be used to differentiate the 
two states. Therefore, HQET has important and unique advantage for this purpose. These 
are the motivation for our approach of using LCQSR in HQET. 

The proper interpolating current J"\>j^^ for the states with the quantum number j, P, 
in HQET was given in [Q. They were proved to satisfy the following conditions 

X y" dtb{x - vt)Jlp^j^{x) (2) 

in the mg oo limit, where ri°'^"'°'3 is the polarization tensor for the spin j state, v 
is the velocity of the heavy quark, g'^^ = g""^ — v'^v^ is the transverse metric tensor, S 
denotes symmetrizing the indices and subtracting the trace terms separately in the sets 
(ai ■ ■ ■ ttj) and {Pi ■ ■ ■ Pj), fpj^ and Hpj^ are a constant and a function of x respectively 
which depend only on P and jg. Because of equations (|lD and (|]), the sum rule in HQET 
for decay widths derived from a correlator containing such currents receive no contribution 
from the unwanted states with the same spin-parity as the states under consideration in 
the rriQ —>■ oo. Starting from the calculations in the leading order, the decay amplitudes 
for finite tuq can be calculated unambiguously order by order in the l/mq expansion in 
HQET. 



2 Sum rules for decay amplitudes 

In the present work we shall confine ourselves to the lowest lying states in the leading 
order of l/mg expansion. Denote the doublet (l+,2+) with ji = 3/2 by (i?i,i?2) and 
the doublet (0"'", 1^) with ji = 1/2 hj {Bq,B[). There are many combinations for the 
coupling constant g-nB^By^ where B^, By belongs to the three doublets (0~,1~), (0^,1'^) 
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and (l+,2+). Due to the heavy quark symmetry and chiral symmetry there exist only six 
independent couphng constants, i.e., three independent couphng constants when B^, By 
belongs to different doublets and another three when B^., By belongs to the same doublet. 
For example, from covariance and conservation of the angular momentum of the light 
component in the mg — *• oo limit, the amplitudes for the decay of Bi, B^ to the ground 
states B, B* in the doublet (0~, 1~) with j£ = 1/2 have the following forms 

M{B, ^ B*n) = I e;vM<lt " ^^?r9t B*) , (3) 

MiB; ^Bn) = I v,uqNt9iB;, B) , (4) 
M{Bl ^ i?*7r) = / le^p^p e*'^v''r]'^'^q?qt,g{B;, B*) , (5) 

where rj^^,, rj^ and are polarization tensors for states 2+, 1+ and 1^ respectively, = 
Ifj, ~ ■ Qv^. I = 1 for charged and neutral pion respectively. As shown in [|1^ the 
constants in (0)-(|^) satisfy 

g{B;, B) = g{B;, B*) = ^ g{Br, B*) . (6) 

We calculate the following six independent coupling constants: gi = gnSiB*, g2 = g-KS'^B, 
gs = gnB[Bi,, gA = gnB*Bi, g5 = gnB[Bi and ge = gjrB*B- Other coupling constants are related 
to them. There are related results for g2, gs and g^ in [6-15]. In gi and g2 is calculated 
with QCD sum rules in the short distance expansion. For the sake of comparison and 
completeness, we also present the sum rules for g2, gs and g^. 

For deriving the sum rules for the coupling constants we consider the correlators 

/ d'x e-''-{n{q)\T (jj^^ ^(0) jj- |0) = (g^^f - ^gfq^'^ I GsM^.^') , (7) 

/ d'x e-^'-%ix{q)\T (jo,+,. (O)^^ . (a:)) |0) = / Gb^bI^,^') , (8) 

j d'x e-'=-(7r(g)|T (v,! (0) ,i (a:)) |0) = q^I Gb',bi,{u^,u;') , (9) 
/ d'x e-^'-{n{q)\T (j{:+,|(0)4t.|(a:)) |0) = 

(gf^r^f - llKgfQt + 9r<l't + \9fqt)) I Gsr^B. {uJ, J) , (10) 

j rf^x e-*^-^(7r(g)|r [jl^.{Q)Jl%dx)^ |0) = le^^'^' q'^v^I Gb'^bA^^uj') , (11) 

j d'x e-^'^-^(7r(g)|r (^^,1 (0)4_^i (x)) |0) = gfJ Gb*b{u:,u:') , (12) 

where k' = k + q, u = 2v ■ k, u' = 2v ■ k' and q^ = 0. The forms of the right hand side of 
(|^-(|1^) are determined by the conservation of angular momentum of the light component 
and the fact that a, P, a and 6 are transverse indices, x — y = vt on the heavy quark 
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propagator. The interpolationg currents are given in M as 



_ . /3 7 5. / 1 



JZ,I = y 4 Kl'H) [V? - -irVt) q , (13) 



where hy is the heavy quark field in HQET and = — 

Let us first consider the function Gb^b* (cj, cu') in (^. As a function of two variables, it 
has the following pole terms from double dispersion relation 

+ -7 + 7T7 ' (17) 



(2A_ 1 -cu')(2A+ 3 -cu) 2A_i-cj' 2A+ s - 



'2 ' ' ' '2 '2 



where /pj^ are constants defined in ([^), Apj^ = mpj^ — uiq. As explained in Section [I|, 
only one state with = 1+ contributes to (|17D as the result of equation This would 
not be true if the last term in ( [T3|) is absent. 

For deriving QCD sum rules we calculate the correlator (^ by the operator expansion on 
the light-cone in HQET to the leading order of 0{l/mQ). The expression for GbiB*{^^,^') 
reads 

^ J^^dt J rfxe-^^-5(-x -t;t)Tr{7*(l + {))75(Di - I7* ^*)(vr(g)|«(0)J(x)|0)} , (18) 

The pion wave function is defined as the matrix elements of nonlocal operators between 
the vacuum and pion state. Up to twist four they are [pi]] : 

< 7r(g)| J(x)7^75w(0)|0 > = -tf^q^ f du e*"«^(^^(M) + x^g^{u) + 0{x^)) 

•J 

+ Ux,-"^) ( due^'^'^^g^^u), (19) 
qx Jo 

< 7r(g)|d(x)i75u(0)|0 > = " / du e'^^^^ifpiu) , (20) 



rriu + nid Jo 

— f TTl^ 

< 7r(g)|d(x)(T^^75n(0)|0 > = i{q^x^ - g^a:^) „, ^ / , / c/m e'"'?'^(^^(u) . (21) 

b{mu + md) Jo 

The wave function (fj^ is associated with the leading twist two operator, gi and g2 corre- 
spond to twist four operators, and ipp and ip^ to twist three ones. Due to the choice of 
the gauge x'^yl^(x) = 0, the path-ordered gauge factor Pexp (igg dux^Af^{ux)) has been 
omitted. 
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Expressing (p!8| ) with the pion hght-cone wave functions, we arrive at 

8 JO JO q ■ V b 

(22) 

where = = 1.76GeV, F-k = ^ = 92MeV for neutral pions. For large euclidean 

values of u and u' this integral is dominated by the region of small t, therefore it can be 
approximated by the first a few terms. 
Similarly, we have: 

Gb'b{oo,u;') = -F^ Hdi /'due^(i-«)1'e^"%'K^p(n) - {q ■ v)[^M + g,{u)]} . (23) 
° 4 JO JO 



^ " 4 Jo JO q ■ V o 

(24) 

Gb*bA^.^') = -i^F^ / dt / due'^^-''^-e'''—u^{^M+t^9i{u)+ 92iu) + -fi^ipM} 

^ o Jo JO q ■ V b 

(25) 

Gb'^bA^^uj') = i^F^ / dt / rfMe^(i-")^e»^M{/x^¥Pp(M) + (g ■ t;)K(M) + tV(w)]} , 
^ 24 Jo JO 

(26) 

Gb*b{^,^') = -F^ / dt / ciMe'(^-")^e'"— {yp,(u) +t2(7i(M) + g2{u) + -/i,yp,(u)} . 

4 Jo JO q ■ V b 

(27) 

After performing Wick rotation and double Borel transformation with the variables uj 
and uj' the single-pole terms in ([T7|) are eliminated. Subtracting the continuum contribution 
which is modeled by the integral in region uj^uj' > uJc, we arrive at: 

9if^if+,l = -^F^^ ^ ~ {uo^Auo)T fo{j^)-^UQgi{uo) + ^Gi{uQ) + ^^^UQ^„{uQ)] , 

(28) 

where uq = j^^jr, T = Ti, T2 are the Borel parameters, Gi{uo) = ug2{u)du 

and /„(x) = 1 — fr- The presence of the factor /„ is the result of subtracting the 

fc=0 

integral s"'e~Tds as a contribution of the continum. In obtaining (^) we have used 
the Borel transformation formula: 3^6°"^ = 6{a — j;). The integral in the function Gi(mo) 
in (^) arises from the factor l/{q ■ v) in (p^). Here we have used integration by parts to 
absorb the factor l/{q ■ v). In this way we arrive at the simple form after double Borel 
transformation. In the following we shall use the same technique to deal with the factor 
{q ■ v)~^ in other sum rules. 
Similarly we have: 

^72/_,l/^.,l = \F^e^^^ {-cp'MT'f li"^) + 2/i.^p(no)T/o(^^) + Ag[iuo)} , (29) 

where ip'^^uo), g[{uo) are the first derivatives of (Pn{u), gi{u) at m = Uq. The appearance of 
derivatives is due to the factor {q- v). 
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93fl 



1 



where GsImq) = Jq° g2{u)du. 



1 



T 



where G^^Uq) = Jq° u'^g2{u)du. 



(30) 



F^e T {uo(Pn{Uo)Tfo{ — ) - -Uo^l(Uo) + -^4(^0) + - finUQifaiUo)} , (31) 



95f+lf+.^ 
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l+A, 

12 Z 



)-4^ (m^i(m)) Ino) 

(32) 



gefl^i = ^F^e T ' {(^^(^o)T/o(^) - ^^i(mo) + ^^3(^0) + ^f^n^aiuo)} ■ (33) 



3 Determination of the parameters 

In order to obtain the couphng constants from (p8|)-(^) we need to use the mass parameters 
A's and the couphng constants /'s of the corresponding interpolating currents as input. 
A_^i/2 and /-,i/2 can be obtained from the results in as A_^i/2 = 0.5 GeV and /-,i/2 — 
0.25 GeV'^/^ at the order as = 0. Notice that the coupling constant /-,i/2 defined in the 
present work is a factor l/\/2 smaller than that defined in flBl. A_|_ 3/2 and A_,_ 1/2 are 
given in and /+,i/2 can be determined from the formulas (34), (27) and (28) of 

reference 0] derived from sum rules for two point correlators. The results are 



Ah 
A. 



,3/2 
,1/2 



0.82 GeV 
1.15 GeV 



f+,3/2 

/+,l/2 = 



= 0.19 ±0.03 GeV^/2 
0.40 ± 0.06 GeV^/2 . 



(34) 



We use the wave functions adopted in [|rT| to compute the coupling constants. Moreover, 
we choose to work at the symmetric point Ti = T2 = 2T, i.e., Mq = | as traditionally done in 
literature |]ll|. We adopt the scale fi = 1.4GeV, at which the values of the various functions 
appearing in ([28|)-(|33|), at uq = |, are: iPn{uo) = 1.22±0.3, (Pp^uq) = 1.142, (Pa-iuo) = 1.463, 
gi{uo) = 3.4 X 10-2 GeV^ Gi{uo) = -4.5 x 10"=^ GeV^, ^3(^0) = -2.0 x 10^^ QgY2^ 
G^iuo) = -1.3 X 10-3 GeV2, g2il) = 0, ^'^{uo) = and g[{uo) = 



4 Numerical results and discussion 

We now turn to the numerical evaluation of the sum rules for the coupling constants. The 
lower limit of T is determined by the requirement that the terms of higher twists in the 
operator expansion is less than one third of the whole sum rule. This leads to T > 1.0 GeV 
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for the sum rules (p8D-(p^). In fact the twist-four terms contribute only a few percent to 
the sum rules for such T values. The upper limit of T is constrained by the requirement 
that the continuum contribution is less than 30%. This corresponds to T < 2.5GeV. With 
the values of pion wave functions at Mq = ^ given above we obtain the left hand side of the 
sum rules (|28| ) - (|33|) as functions of T. The results are shown in FIG. 1. Stability develops 
for the sum rule (p8D-(p3D in the region 1.0 GeV <T<2.5 GeV. Numerically we have: 



9if-,hf 

92 f- 

93 f I 
9.fl 

95 f+ 



= -(0.16 ±0.02 ±0.03) GeV^ 
i/^ 1 = (0.36 ±0.02 ±0.03) GeV^ 
1 = (0.13 ±0.01 ±0.02) GeV^ 

2 

3 = (0.10 ±0.01 ±0.02) GeV^ 

2 

i/+ 3 = (0.20 ±0.01 ±0.02) GeV^ 



96fl,i = (0.20 ± 0.01 ± 0.02) GeV^ 



(35) 
(36) 
(37) 
(38) 
(39) 
(40) 



where the first error refers to the variations with T in this region and the second error 
takes into account the uncertaity in Uc- And the central value corresponds to T = l.SGeV 
and Uc = 3.0GeV for (||)-(|3|). For the sum rule (H) we use = 2.4GeV. 

With the central values of f 's in (B^ we get the absolute value of the coupling constants: 



91 

92 
93 
94 
95 
96 



-(3.6 ±0.4 ±0.6) 

3.6 ±0.3 ±0.6 , 
0.83 ±0.1 ±0.2 

2.7 ±0.2 ±0.5 
2.6 ±0.2 ±0.4 
3.1 ±0.3 ±0.6 



GeV 



GeV^^ 
GeV-3 
GeV-i 
GeV-i 



(41) 
(42) 
(43) 
(44) 
(45) 
(46) 



We have not included the uncertainties due to f's here 
and (|30D -(|33 D are also sensitive to the value of (Pn{^)- 



Note that the sum rules 
For example, gi = —4.83 with 



1.71 which is one of the values used in 
The couplings of pions with heavy mesons obey symmetry relations in the leading order 
of HQET. Only in HQET the ambiguity due to the presence of two distinct 1+ states can 
be resolved. We use light-cone QCD sum rules to calculate all the six independent strong 
coupling constants gi-g^ for the lowest three doublets. The sum rules in HQET for g^ and 
gQ are consistent with those for g and g' in |]1T], |13| in the limit of tuq oo with gs = ^ 



and ge 



2g. 



The sum rule 



for g2 can be compared to that for h in [13| with g2 



2h 



from definition. In terms of h our value is 0.25, which is only one half of that given in |T3 
Our result gi = 3.6 is slightly smaller than that from the short-distance QCD sum rule 
The calculation for g^ and g^ is new. 

With the CLEO measurement of B{D*^ —>■ D^'-f) ||18[, a model independent extraction 
of g has been performed recently [0. Two possible solutions of g { g 
notation ) and /3 are found from fitting to the experimental data, either g 
(3 = 0.85l°:?t°:?GeV-i or g = 0.761°:°^°;?, (3 = 4.90l°j;^:?GeV-i. The g 
is excluded by the experimental limit r{D*^) < 0.13MeV The coupling g from the 



= ^ge m our 

n O7+0.04+0.03 
U-^' -0.02-0.02) 

: 0.76 solution 
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-0.40 R,|T21P 



LCQSR approach is (7 ~ 0.21 
favored solution g = 0.27. In a recent work 
lowest excited heavy meson doublets (2+, 1 
the leading order coupling gi was calculated. The resulting pionic decay widths of 1^2(2460) 
and Di(2420) are in agreement with the experimental data p2| . 



which is consistent with the experimentally 
1| the mixing effect of the two 1+ states in the 
and (1"'",0~*") and the O^I/tjiq) correction to 
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Figure Captions 



FIG. 1. Dependence of fifjgi-e in (|2^)-(|3^ on the Borel parameter T for 
the continuum threshold Wc = S.OGeV. The dotted, solid, short-dashed, 
long-dashed, dot-dashed and intermediate-dashed curves correspond to 
the coupling constant gi-ge respectively. 



9 



